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INTRODUCTION 


The  success  of  a  farm  enterprise  in  arid  regions  often  depends 
on  the  ability  of  the  farmer  to  adapt  dependable  and  profitable  farm- 
ing practices.  Within  such  regions  one  of  the  primary  climatic  ele- 
ments influencing  crop  production  is  precipitation.  The  objective  of 
this  study  was  to  investigate  a  method  to  characterize  annual  precipi- 
tation and  crop  yields. 

Evaluation  of  the  dependability  or  conversely  the  risk  associated 
with  a  crop  at  a  given  location  in  the  past  has  often  followed  regression 
techniques.  An  example  of  this  approach  is  found  in  Mathews  and  Brown 
(7)«3  In  this  instance,  the  procedure  consisted  of  deriving  a  linear 
relationship  between  yields  of  crops  and  variables  such  as  soil  mois- 
ture or  precipitation.  An  estimate  of  the  frequency  of  occurrence  of 
various  yield  levels  for  an  area  was  then  determined  by  applying  the 
derived  relationship  to  point  rainfall  data  from  that  area.  Charac- 
terization of  annual  precipitation  amounts  has  proceeded  little  past 
the  determination  of  means  and  coefficients  of  variations.  However, 
Potter  (10)  has  shown  that  the  extreme  value  distribution  may  be  used 
to  characterize  the  occurrence  of  annual  precipitation  amounts. 

One  of  the  limitations  of  the  regression-approach  to  determine 
the  dependability  of  crop  production  lies  in  the  rationale  for  the  re- 
lationship between  yields  and  soil  moisture  and/or  rainfall.  Yields 
are  influenced  by  phenomena  not  necessarily  related  to  soil  moisture 
or  precipitation  conditions,  e.g.,  animal  or  insect  infestations,  wind, 
or  hail.  Further,  the  utilization  of  this  approach  requires  laborious 
tabulations  and  computations. 
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With  the  above  considerations  in  mind,  an  improvement  in  tech- 
nique for  evaluating  or  characterizing  the  dependability  of  a  crop 
or  cropping  practice  or  precipitation  should  possibly  employ  a  proce- 
dure having  the  following  properties : 

i)  Have  a  rational  model, 
ii)   Give  consideration  of  all  factors  which  influence  the 
phenomenon . 
iii)   Provide  a  method  which  is  simple  computationally. 

During  recent  years  the  statistical  distribution  of  the  extreme 
values  in  a  sample  has  been  utilized  with  success  in  many  diverse  ap- 
plications (5).  As  applied  to  the  largest  extreme,  the  fundamental 
theorem  has  been  stated  as  follows  (h) : 

"In  a  set  of  N  independent  extremes  y]_,  y2,  ,  y^,  each 

being  the  extreme  of  one  of  N  sets  of  n  observations  each  of 
an  unlimited,  exponentially-distributed  variable,  as  both  N 
and  n  grow  large  the  cumulative  probability  that  any  one  of 
these  N  extremes  will  be  less  than  any  chosen  quantity,  y, 
approaches  the  double  exponential  expression 


_  a  (y-y) 
e 
F(y)  =  e"  [ij 

where  y  is  the  mode  of  the  set." 

From  fjQ ,   the  average  return  period  of  a  quantity  y  being  equalled  or 
exceeded  is  Tm  =  l/  (l-F(y)  ).  fol 

A  crop  yield  in  any  one  year  at  a  particular  location  is  an 
observable  quantity  and  its  magnitude  is  governed  by  the  interaction 
of  environmental  components.  Such  components  include  soil  properties, 
farming  practices,  precipitation,  temperature,  insect  infestations 
and  other  factors.  Conceptually  then,  the  observed  yield  comes  from 
a  large  population  of  possible  values  ranging  from  zero  to  some  upper 
bound.  It  seems  reasonable  that  the  conceptual  population  of  yields 
should  be  distributed  approximately  in  an  exponential  manner,  being 
bounded  from  below  and  above.  Under  dryland  farming  conditions  most 
observed  annual  yields  are  far  removed  from  the  upper  bound  and,  to 
a  lesser  extent,  from  the  lower  bound  of  zero.  Recorded  zero  yields 
may  be  indicative  of  a  crop  not  warranting  harvest  and  in  this  event 
would  not  strictly  be  a  zero  yield.  Thus,  with  reservations  the  con- 
ceptual population  of  annual  yields  may  approximate  the  unlimited  ex- 
ponential distribution  required  by  the  fundamental  theorem.  The 
observed  yield  for  a  given  location  and  year  can  be  treated  as  a  maximum. 
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The  same  reasoning  applies  to  observed  annual  precipitation. 
The  observed  annual  amount  is  from  a  large  conceptual  population  and 
can  be  treated  as  a  maximum.  Here  again  the  population  may  be  bound- 
ed by  zero  and  some  unknown  upper  bound.  However,  the  lower  bound, 
for  the  observational  period  considered  is  rarely  if  ever  attained. 
It  is  not  known  whether  the  upper  bound  is  ever  attained.  Again  it 
is  hypothesized  that  the  fundamental  theorem  is  satisfied,  possibly 
to  a  closer  extent  than  with  annual  crop  yields. 

This  study  is  an  empirical  investigation  pertaining  to  the  val- 
idity of  applying  extreme  value  theory  to  annual  yields  and  precipi- 
tation amounts,  i.  e.,  whether  the  model  is  reasonable.  Data  from 
several  locations  were  utilized  for  this  purpose  (l,  2,  6,  8,  9). 
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APPLICATION  OF  THEORY 


GENERAL  RELATIONSHIPS 


Detailed  information  on  fitting  equation  /!/  to  observed  data 
can  "be  found  in  (3,  ^,  5).  Equation  QJ  can  be  linearized  by  taking 
an  iterated  natural  logarithm  of  both  sides,    i.e., 

z   =   -In  (-In  F  (y)      )      -   a  (y  -  y)  ffl 

Introducing  two  theoretical  quantities,  see  (K) } 

a  =  -^Nj_  and  y  =  y  -  Sv  J^_     , 
by  y       -N 


equation  [h]   becomes 


y   =   y  +   S           * ±» 

Dividing  J^J  by  y,   one  has 

(  y  )         1  ,  (sy  )       (  z 

-  yN   )    , 

y  y       ^N 

where 

(   y  )      coefficient  of  variation  of 


5J 


3J 


y 


the  series  of  extremes,  C.  V. 


z   is  defined  by  /3/,°NN  and  y^  are   tabulated  quantities   (k) ,see 
Table  2. 


PLOTTING,  COMPUTATIONAL,  AND  FITTING  PROCEDURES 

For  an  examination  of  how  the  theory  represents  the  observed 
series,  it  is  convenient  to  plot  the  data  and  to  fit  equation  [6]  to 
the  plotted  points.  The  method  of  procedure  is  as  follows: 

(a)  The  set  of  ^  _   values  are  ordered  according  to  magnitude 

with  the  largest  value  ^  first,  the  next  largest  value  second, , 

the  smallest  value  in  the  N   position. 


(b)  Assign  a  rank,  m,  to  the  above  ordered  series,   i.  e.,  the 
first  value  has  a  rank  of  one,  the  next  largest  value  a  rank  of  two, 
,  and  the  smallest  value  a  rank  of  N. 


(c)  For  each  ordered  value  calculate  F  ( = — )  as  follows: 

where  m  is  the  rank  of  a  particular  ( — == — )  value  and  N  is  the  total 
of  observations. 

If  special  probability  paper  is  available  (5)  the  points 
(  y/y?  F  (y/y)  )  are  plotted.  However,  if  such  paper  is  not  avail- 
able it  can  be  constructed  or  values  of  (y/y>  -In  (-In  F  (y/y)  )  )  may 
be  plotted  on  linear  coordinate  paper. 

Calculation  of  the  quantity  (Q-y/—)  is  carried  out  with  either  the 
ranked  or  unranked  data.  Thus, 

(Sy/y)2  ■   —?_       S  (y/y)2    "    _N_ 
J  N  -  1  N  -  1 


and  C.  V.  =  Sy/_  =   /  _1_       (  S  (y/y)2   -  N  V 


For  fitting  equation  /6j  to  the  data  it  is  convenient  to  utilize 
determined  values  of  y  and  <^N,  the  reduced  means  and  standard  devia- 
tions, respectively,  for  obtaining  a  "least  squares"  type  fit.  Values 
associated  with  N  =  20  to  N  =  59 >  inclusive,  are  given  in  Table  2.  Two 
points  of  equation  /6/,  utilizing  Table  2,  can  be  located  on  the  graph 
of  (y/y),  F  (y/y)  for  orienting  the  straight  line. 


APPLICATION  OF  PLOTTING  AND  FITTING  PROCEDURE 

A  detailed  example  of  the  plotting  and  fitting  procedures  will 
be  of  value  for  possible  applications.  For  this  purpose  V7,  (N  =  ^l) , 
years  of  annual  precipitation  data  (Jan.  1  -  Dec.  31)  from  Akron, 
Colo.,  were  used,  see  Table  1.  Values  shown  in  Table  1  for  each  of 
the  ranked  precipitation  amounts  were  obtained  as  follows: 

N  =  ^7 

P  =  2]  P  /  N  =  8°3-6      =   17.097  ins. 


N 

_ 

803.6 

W 

N  + 

1 

-  m 

F  (_£j   =   N+  1  -  m    =   hQ   -  m    ,      ..  .  ,       r-j 
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Calculation  of  Sp/_  ,    coefficient  of  variation,   C.  V.,  is  found  as 
follows:         P 

2    (P/p)2       =       li9.35786l. 

(Sp/-)2    =       ^9.357861         -        J£L 
*  p  1+6  1+6 

(Sp/p)        =       C.   V.        =        0.2261+ 

Thus  the   fitted  equation  is 

(   P  )       =       1  +  0.2261+       (z   -  yk  )  /|7 

P  ^N 

With  N  =  1+7,  values  of  yN  and  <=^N  from  Table  2  are  found  to  be 
0.5^73  and-  1.1557,  respectively.  Two  points  on  equation  [Gj ,  chosen 
arbitrarily,  can  now  be  found  as  follows: 

Assume  that  F  (p/p),   =   0.990, 

then  zx,        =        -In  (-In  (0.990)      )      =     1+.60517 

and  assume  that  F   (P/p)2     =     0.800, 

then  z       =      -In     (-In  (0.800)      =     I.5OO58. 

Substituting  the  values  of  ^N,      and  <t~n,   and  the  calculated  z 
and  z0  values   into  equation  [6j  gives  the   following: 

P_\      =       1  +   0.2261+       (k. 60517   -  0.5^73) 


(*) 


p;  1.1557 

p_\    =  1.79^9 
p 


and 


.♦). 


1.  +    0.2261+        (I.50058   -  0.5^73) 

1.1557 


1873 

The  "least  squares"  equation  will  pass  through  the  points 


^)i     FWj    md        ((p/?)2-    F(- 

i.  e.,  (1.79^9,  O.990)  and  (1.1873,  0.800). 
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TABLE  1.--AKNUAL  PRECIPITATION  DATA,  AKRON,  COLO. 


p 

ins. 

P 
P 

Rank 

m 

F(!) 

p 

P 
ins. 

p 

p 

m 

F  (   I) 

P 

26.5 

1.550 

1 

0.979 

16.7 

0.976 

25 

O.479 

25.0 

1.1*62 

2 

0.958 

16.7 

0.976 

26 

O.458 

23.1 

1.351 

3 

0.937 

16.6 

0.971 

27 

0.437 

22.5 

1.316 

4 

0.916 

15.7 

0.918 

28 

0.4l6 

22.3 

1.301+ 

5 

0.895 

15.6 

0.912 

29 

0.395 

22.1 

1.293 

6 

0.875 

15.6 

0.912 

30 

0.374 

21.6 

1.263 

7 

0.854 

15.5 

0.906 

31 

0.353 

21.  4 

1.251 

8 

0.833 

15.5 

0.906 

32 

0.332 

20.7 

1.211 

9 

0.812 

15.1+ 

0.900 

33 

0.311 

20.  4 

1.193 

10 

0.791 

15.1 

0.883 

3^ 

O.290 

20.2 

1.181 

11 

0.771 

14. 5 

0.81+8 

35 

O.269 

19.6 

1.146 

12 

0.750 

14.1+ 

0.81+2 

36 

0.248 

19.1 

1.117 

13 

0.729 

li+.i 

0.825 

37 

O.227 

19.0 

1.111 

14 

0.708 

13.7 

0.801 

38 

0.206 

18.9 

1.105 

15 

0.687 

13.!+ 

O.78I+ 

39 

O.185 

18.7 

1.094 

16 

0.666 

12.7 

O.7I+3 

1+0 

0.164 

18.5 

1.082 

17 

0.61+6 

12.6 

0.737 

1+1 

O.143 

18.2 

1.061+ 

18 

O.625 

11.8 

O.690 

1+2 

0.122 

18.2 

1.061* 

19 

0.604 

11.7 

0.681+ 

1+3 

0.101 

18.1 

I.058 

20 

O.583 

11.1+ 

0.666 

1+1+ 

0.080 

17.7 

1.035 

21 

O.562 

11.0 

O.6I+3 

^5 

0.059 

17.5 

I.023 

22 

0.51+2 

10.7 

0.626 

1+6 

0.038 

17.4 

I.017 

23 

0.521 

9.7 

O.567 

47 

0.017 

16.8 

O.982 

21* 

0.500 

TABLE  2  -REDUCED  MEANS  AND  STANDARD  DEVIATION1 


N 

y 

N 

<r>  u 

N 

^N 

cr^N 

20 

0.5236 

1.0628 

ho 

O.5I+36 

1.11*13 

21 

0.5252 

1.0696 

1*1 

O.5I+1+2 

1.11*36 

22 

0.5268 

1.0751* 

1*2 

O.5I+1+8 

1.11*58 

23 

0.5283 

1.0811 

h3 

0.5^53 

1.11*80 

2k 

0.5296 

1.0861* 

1+1+ 

O.5I+58 

1.11+99 

25 

0.5309 

I.0915 

^5 

0.51*63 

1.1519 

26 

0.5320 

I.0961 

1+6 

O.5I+58 

1.1538 

27 

0.5332 

1.1001+ 

*7 

0.5^73 

1.1557 

28 

o.53^3 

1.101*7 

1+8 

0.5^77 

1.1571* 

29 

0.5353 

1.1086 

1+9 

0.51*81 

1.1590 

30 

0.5362 

1.1124 

50 

0.51*85 

1.1607 

31 

0.5371 

1.H59 

51 

O.5I+89 

1.1623 

32 

0.5380 

1.H93 

52 

0.51*93 

1.1638 

33 

O.5388 

1.1226 

53 

0.51*97 

1.1658 

3^ 

0.5396 

1.1255 

5^ 

O.5501 

1.1667 

35 

0.51*03 

1.1285 

55 

O.550I+ 

1.1681 

36 

0.51*10 

1.1313 

56 

O.5508 

1.1696 

37 

0.51+18 

1.1339 

57 

0.5511 

1.1708 

38 

0.5if2i+ 

1.1363 

58 

0.5515 

1.1721 

39 

0.51+30 

1.1388 

59 

0.5518 

1.1731* 

*  Table  from  p.  65  (4)      Used  with  the  permission  of  the  authors  and  publishers. 
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AGREEMENT   BETWEEN  THEORY  AND  OBSERVATIONS 


Figure  1  shows  the  accumulated  probability  F  (  /p)  and  F  ("Vy) 
corresponding  to  levels  of  (  /p)  and  (  /y)  at  selected  locations. 
Shown  in  this  figure  are  the  plottings  of  both  precipitation  and  crop 
data  from  6  locations  for  given  lengths  of  record,  N.  Average  pre- 
cipitation, P,  is  in  inches  and  yields  are  in  bushels  per  acre  ex- 
cept cotton  which  is  in  pounds  per  acre.  Coefficients  of  variation, 
C .  V. ,  are  also  given  for  each  series  of  data.  The  broken  lines  on 
each  plot  indicate  the  "50  per  cent"  confidence  bands.  Procedures  for 
construction  of  these  bands  can  be  found  in  (3,   ^ }    5). 

It  has  been  stated  (k)  that  if  two-thirds  of  the  observations, 
as  plotted  on  extreme  value  probability  paper,-  fall  within  the  confi- 
dence band,  the  observed  series  may  be  considered  to  be  represented 
adequately  by  the  theory  of  extreme  values.  For  the  situations  stud- 
ied here  this  criterion  is  satisfied.  As  would  be  deduced  from  the 
fundamental  theorem,  the  series  of  continuous  wheat  yields  in  which 
zero  values  are  encountered  indicated  the  poorest  fit.  Conversely, 
for  situations  in  which  zero  values  seldom,  if  ever,  occur,  i.e.,  fal- 
low wheat  yields  and  annual  precipitation,  a  better  fit  was  evident. 


DISCUSSION 


UTILITY  OF  THE  EXTREME  VALUE  THEORY 

The  simplicity  of  the  present  technique  is  evident  in  that  only 
the  quantity  Sy  and  plotting  positions  need  to  be  calculated  from  the 
data.  y  Construction  of  special  probability  paper  makes  the 
plotting  procedure  simple. 

Having  a  description  of  the  past  population  of  crop  yields  or 
precipitation  amounts,  it  would  be  possible,  if  the  observed  popula- 
tion is  characteristic  of  the  larger  parent  population,  to  predict 
selected  crop  or  precipitation  characteristics  with  an  associated  prob- 
ability statement.  It  should  be  pointed  out  that  the  "classical"  ap- 
proach such  as  deriving  a  relation  between  yields  and  precipitation  in 
order  to  predict  future  yields  must  utilize  some  prediction  or  proba- 
bility description  of  future  precipitation.  With  the  present  approach 
the  intermediate  regression  step  does  not  exist  as  an  estimate  can  be 
made  directly  from  the  distribution  of  past  yields  wherein  all  factors 
influencing  yield  have  manifested  themselves. 

An  examination  of  January  to  December  precipitation  and  July  to 
June  precipitation  led  to  the  hypothesis  that  any  annual  total  would 
yield  similar  results.  Further,  monthly  and  seasonal  amounts  of  pre- 
cipitation were  found  amenable  to  description  by  the  technique. 
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Figure  1.    Relationship  between  accumulative  probability,  F,  and  precipitation  and  yield  ratios  of  annual  values 
to  long  time  means. 
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AN   EXAMPLE  OF   CROP   DEPENDABILITY  AS  INFLUENCED 
BY   FALLOW  OR  CONTINUOUS  CROPPING: 

The  dependability  of  a  crop  can  be  measured  by  the  frequency  of 
occurrence  of  yields,  for  a  selected  length  of  record,  within  arbitrary 
intervals,  e.g.,  Y  <  5  bu./ac,  5  bu./ac.^Y  <15  bu./ac.  From  the 
development  presented  such  determinations  are  readily  obtained  by  the 
following  expression, 

Pr   |  L  ^  Y  <  U  I  =      Pr  |  Y  <  U  |    -    Pr   I  Y  <  L 


where  Pr  is  the  probability  of  the  yield  level  Y  falling  between  arbi- 
trary levels  L  and  U.  Each  of  the  expressions  on  the  right-hand  side 
are  directly  read  from  the  pertinent  graph  of  Figure  1. 

As  an  example  of  evaluating  crop  dependability,  winter  wheat 
yields  at  Akron,  Colo. ,  and  Woodward,  Okla.,  are  contrasted  with  respect 
to  a  continuous  and  summer  fallow  cropping  practice. 

In  order  to  evaluate  yield  interval  probability  statements  uti- 
lizing Figure  1  the  interval  limits  will  have  to  be  changed  to  ones  for 
the  ratio  Y/y»  F°r  example  at  Akron,  the  average  Y,  for  continuous 
wheat  is  7*532  bu./ac.  Thus  values  of  Y  equal  to  5  bu./ac.  and  15  bu. 
/ac.  correspond  to  values  of  Y/y  of  0.66^  and  1.991>  respectively. 
Probability  values  are  derived  as  follows: 


Pr    J  Y  <   5  j  =     Pr    J  Y/y  <  0.66U  j    =  O.kk 
Pr   I  5  <  Y  <  15  I  —  Pr   j  0.66k   <    Y/v  4  1-99l 


Pr    j  y/y  <  I.991  I  -   Pr   j  Y/y<  0.66k   j  =  0.81  -  0.^=0.37 


An  examination  of  Table  3  clearly  shows  why  fallow  cropping  is 
preferred  at  Akron  and  continuous  cropping  is  preferred  at  Woodward. 
It  must  be  remembered  that  fallow  cropping  cuts  the  acreage  harvested 
annually  in  half.  From  the  first  line  of  Table  3  it  will  be  noted  that 
the  probability  of  obtaining  a  yield  less  than  5  bu./ac.  at  Akron  for 
continuous  cropping  is  O.kk.  With  fallow  wheat  it  is  0.09.  The  hazard 
of  crop  failure,  therefore,  is  almost  5  times  greater  with  continuous 
wheat.  By  contrast,  at  Woodward,  the  probability  of  obtaining  a  crop 
failure  is  at  the  low  level  of  0.02  for  both  methods  of  cropping. 
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TABLE  3   -PROPORTION    OF   THE   TIME   WHEAT 
YIELDS  FALL  IN  VARIOUS  YIELD  INTERVALS 


AKRON, 

COLO. 

WOODWARD,  OKLA. 

YIELD  INTERVAL 

Continuous 

Fallow 

Continuous 

Fallow 

Bu.   per  acre 

<     5 

O.i+U 

0.09 

0.02 

0.02 

^    5 

<15 

0.37 

0.35 

0.1+1 

0.30 

^>  15 
«^25 

o.ih 

0.31 

0.1+2 

O.38 

^25 
<  35 

0.03 

0.16 

0.12 

0.19 

^35 

0.02 

0.09 

0.03 

0.11 
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SUMMARY  AND  CONCLUSIONS 

For  the  data  considered  in  this  study  the  extreme  value  theory 
satisfactorily  represents  the  distribution  of  annual  crop  yields  and 
annual  precipitation  amounts.  The  computational  and  graphical  proce- 
dure is  simple.  It  is  possible  to  characterize  selected  crop  or  pre- 
cipitation measures  by  their  probability  of  occurrence. 

The  extreme  value  technique  presented  herein  appears  to  offer  a 
technique  adapted  for  an  extensive  survey  of  cropping  risks,  as  well 
as  expected  precipitation  amounts.  A  study  of  this  type  would  be  use- 
ful in  depicting  the  hazards  of  crop  production  in  areas  such  as  the 
Great  Plains.  For  some  purposes  there  may  be  advantages  in  working 
directly  with  the  distribution  of  crop  yields  and  allowing  all  factors 
which  influence  them  to  be  taken  into  account. 

The  marked  difference  between  the  annual  precipitation  and  yield 
graphs  indicates  that  annual  precipitation  for  the  situations  studied 
here  Is  not  a  satisfactory  indicator  of  the  variability  of  crop  yields. 
Based  on  the  levels  of  the  coefficient  of  variation,  annual  precipita- 
tion showed  less  variability  than  annual  yields. 
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